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^ ! ABSTRACT 

Aims. The homogeneous spectroscopic determination of the stellar parameters is a mandatory step for transit detections from space. 
Knowledge of which population the planet hosting stars belong to places constraints on the formation and evolution of exoplanetary 
systems. 

Methods. We used the FLAMES/GIRAFFE multi-fiber instrument at ESO to spectroscopically observe samples of stars in three 
Ca/?o7YExoplanet fields, namely the LRaOl, LRcOl, and SRcOl fields, and characterize their stellar populations. We present accurate 
atmospheric parameters, T e fr, log g, [M/H], and [o/Fe] derived for 1 227 stars in these fields using the MATISSE algorithm. The latter 
is based on the spectral synthesis methodology and automatically provides stellar parameters for large samples of observed spectra. 
We trained and applied this algorithm to FLAMES observations covering the Mg i b spectral range. It was calibrated on reference 
stars and tested on spectroscopic samples from other studies in the literature. The barycentric radial velocities and an estimate of the 
CZ) , V sin i values were measured using cross-correlation techniques. 

Results. We corrected our samples in the LRcOl and LRaOl CoRoT fields for selection effects to characterize their FGK dwarf 
stars population, and compiled the first unbiased reference sample for the in-depth study of planet metallicity relationship in these 
CoRoT fields. We conclude that the FGK dwarf population in these fields mainly exhibit solar metallicity. We show that for transiting 
planet finding missions, the probability of finding planets as a function of metallicity could explain the number of planets found in 
the LRaOl and LRcOl CoRoT fields. This study demonstrates the potential of multi-fiber observations combined with an automated 
■ classifier such as MATISSE for massive spectral classification. 

C^") ' Key words. Techniques: spectroscopic - Stars: fundamental parameters, planetary systems, general 
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\ 1 . Introduction photometric precision for more than 120000 stars. However, an 

, , , , — — — - — I accurate knowledge of the fundamental parameters of the stars 

Since, the discovery of the first exoplanet by | Ma V or & Quelozj obseryed by ^ satemte h mandatory to fully exploit Ms pho . 



X 



(119951) . large-scale spectroscopic surveys for finding planets tometric database 
have gathered thousands of high resolution spectra. Atmospheric 

parameters and chemical compo s ition for these sampl e s wer e We performed a first-order spectral classification of the stars 



Vh 1 determined by IValenti & Fischer! d2005l) or lSousa et all d2008). in some of the CoRo r/Exoplanet fields using broad-band pho- 
I They explored the links between the hosting stars and field tometry obtained at the INT/La Palma. The coordinates and 
stars parameters. This requires a good characterization of the magnitudes as well as these photometric spectral types and lu- 
fields in which planets are s earched for. The CoR oT (Convection minosity classes of the CoRoT stars are available from Exo-Dat 
Rotation and Transits, see Bagl in et al.l 120061 for full details) dDeleuil et al.||2009t |Meunier et al.||2007|) and used for the target 
space mission has obtained light curves with very high relative selection and the precise placement of its photometric masks. 
However, this spectral classification presents some uncertainties, 



Send offprint requests to: Jean-Christophe Gazzano, e-mail: jean- e.g., the unknown star's reddening, chemical abundances, and 

christophe.gazzano@oamp.fr potential binarity. Combining this photometry with intermediate 

* Based on observations collected with the GIRAFFE and resolution spectroscopy would help in the determination of the 

UVES/FLAMES spectrographs at the VLT/UT2 Kueyen telescope physical parameters of the stars 
(Paranal observatory, ESO, Chile: programs 074.C-0633A & 081.C- 

0413A). In this context, we present a fully automated and homoge- 

** Full tables|3[9HH and[B]are only available in electronic form at neous determination of atmospheric parameters from interme- 

the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via diate resolution spectra obtained with the FLAMES/GIRAFFE 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-t-A/523/A91 multi-object facility. This work is the second scientific objective 
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of the study presented by lLoeillet et all (l2008h . hereafter iLfjH 
In that paper, the authors demonstrated the capability of multi- 
fiber instruments to find planetary candidates by radial velocity 
techniques. We extended this study by observing a new sample 
of stars with the same configuration in May-June 2008, to per- 
form a homogeneous spectral characterization of the fields ob- 
served by CoRoT. For that purpose, we adapted an algorithm 
originally developed for the spectroscopic analysis of data to 
be obtained with the Radial Velocity Spectrometer of the Gaia 
mission: MATrix Inversion for Spectral Synth Esis, MATISSE 
dRecio-Blanco et all l2006t iBijaoui et all 120081) . We measured 
the effective temperature (T e ff), the surface gravity (log g), the 
overall metallicity ([M/H]), and the a-enhancement ([ayFe]) of 
1 227 stars, the barycentric radial velocity of 1 534 stars, and es- 
timated the projected rotational velocity ( V sin z) of 1 604 stars 
located in the LRaOl, LRcOl, and SRcOl fields. The whole data 
set, observational set-ut and date, and the derived physical pa- 
rameters are available to the community through the CoRoT 
database, Exo-DaQ We applied the relation linking the prob- 
ability of finding plan ets with the metallicity of the host stars 
dUdrv & Santosll2007l) to the de-biased sample. 

The structure of this paper is as follows. Section|2]describes 
the observations, the instrument setup, and the target selection. 
Section[3]continues with data reduction and processing. The au- 
tomatic algorithm, its implementation, and limitations are devel- 
oped in Sect. [4] The results and a discussion can be found in 
Sects. [5]and[6j respectively. 

2. Observations 

In January 2005 and May-June 2008, we obtained 13 half-nights 
in visitor mode to perform spectroscopic observations with 
the FLAMES multi-object facility coupled with the GIRAFFE 
and UVES spectrographs (programs 074.C-0633A and 081.C- 
0413A). The instrument is mounted on the 8.2 m Kueyen tele- 
scope (UT2) based at the ESO-VLT. 

The FLAMES/GIRAFFE observations were performed in the 
MEDUSA configuration, using the HR9B spectral domain. This 
setup covers about 200 A centered at 5258 A around the Mg i b 
lines, with an intermediate resolving power (R = 25 900) and 
a CCD pixel sampling of 0.05 A. This instrumental configura- 
tion was selected in 2005 for it contains man y thin spectral line s 
leading to a good radial velocity accuracy (Ro ver et alj |2002). 
This wavelength range is also interesting for the determination 
of spectroscopic parameters since it contains many metallic lines 
that can constrain the temperature and metallicity, and some ion- 
ized spectral lines that can help to constrain the surface grav- 
ity. As a part of the radial velocity follow-up of CoRoT exo- 
planet candidates, we observed 7 additional stars, at a higher 
spectral resolution and across a wider wavelength range with the 
UVES/FLAMES facility. We used the red arm of the spectrograph 
at a central wavelength of 5800 A, covering about 2000 A with 
a resolving power of about 47 000. 

The 2005 campaign was dedicated to the radial velocity 
follow-up of selected CoRoT stars in the so-called anticenter 
direction field (LRaOl: Long Run Anticenter 01), at Galactic 
coordinates I ^ 212.2°, b ^ -1.9°. Full description of the tar- 
get selection and observational strategy for this first campaign 
can be found in Sect. 2 of lL08l The targets observed during 
the 2008 campaign belong to the SRcOl (Short Run Center 
01; I - 36.8°, b = -1.2°) and LRcOl (Long Run Center 01; 
I 37.7°, b ^ -7.5°) CoRoT fields, close to the Galactic center 

1 http://lamwws.oamp.fr/exodat/ 



Table 1. Priority Criteria. 9 is the highest priority, 1 the lowest.. 



Priority 


Sp. Types 


Lum. Class 


Magnitude 


Variability 


9 


F,G,K,M 


IV, V 


r' < 14 


no 


6 


F,G,K,M 


IV, V 


14 < r' < 15 


no 


3 


F,G,K,M 


rv,v 


14 < r' < 15 


yes/no 


1 


All Other Stars 


r' > 15 


yes/no 



Notes. This selection was performed using the information available in 
Exo-Dat. The last column is a fl ag on the probability for the star to be 
variable (Debosscher et al. 2007). 



direction. The observation strategy for the 2008 campaign was 
slightly different from the 2005 one since the main purpose of 
the GIRAFFE observations was to characterize the stellar popu- 
lation observed by CoRoT, which was only a secondary objective 
for the 2005 campaign. 

For the 2008 campaign, each observed field was centered 
around the CoRoT planetary candidates observed with UVES. 
The FLAMES/GIRAFFE fiber allocation was performed taking 
into account the instrumental constraints (magnitude limits, fiber 
positions, etc.) and according to four different levels of priorities, 
as described in Table Q] Using the spectral classification avail- 
able in Exo-Dat, we gave the highest priority to solar-type dwarf 
and sub-giant stars since these were the stellar populations we 
wished to characterize. High priorities were assigned to bright 
targets to ensure good signal-to-noise ratios (SNRs) in short ex- 
posure times. About 4-5 fibers were used to register the sky sig- 
nal. These selection criteria were chosen to be very similar to 
the 2005 ones to ensure the homogeneity of the whole sample. 
Bad weather conditions and instrumental issues affected the sec- 
ond observing campaign. Out of the 20 planned configurations, 
we managed to observe only 10 FLAMES fields, as listed in the 
journal of the observations reported in Table [2] 

In total, we obtained 1241 spectra of stars located in the 
SRcOl and LRcOl CoRoT fields. These new spectra were ana- 
lyzed together with the 772 spectra of LRaOl stars. Hereafter, 
we call the pointing direction related to the LRcOl and SRcOl 
fields center and related to the LRaOl field anticenter. 

Table 2. Journal of the 2008 observations with the coordinates 
of the field centers and the adopted exposure time. 



Field Obs. Date R.A. Dec Exp.T. N 







(hh:mm:ss) 


(dd:mm:ss) 


(Seconds) 




LRc01.12 


2008-06-09 


19:30:14 


-00:00:10 


4290 


128 


LRcOl. 02 


2008-06-09 


19:23:34 


+01:16:36 


2700 


129 


LRcOl.04 


2008-06-09 


19:25:36 


+01:26:46 


3600 


128 


LRc01-05 


2008-05-28 


19:24:19 


+00:49:41 


2401 


128 


LRcOl. 06 


2008-06-09 


19:24:17 


+00:47:09 


3600 


126 


LRc01-09 


2008-06-12 


19:26:57 


+00:38:13 


2400 


129 


SRcOl.01 


2008-06-11 


19:04:14 


+02:02:11 


2700 


128 


SRcOl. 03 


2008-06-12 


19:03:49 


+03:24:50 


3600 


129 


SRc01-05 


2008-06-11 


19:03:00 


+02:47:01 


2700 


129 


SRc01.07 


2008-06-11 


19:03:43 


+02:58:54 


2700 


128 



Notes. The last column gives the nu mber of spectra registered per field. 
For the 2005 observations see lL08l 



3. Data processing 

The 2005 frames were reduced using the G IRAFFE BaseLine 
Data Reduction Software (girbldrs vl.12, see lRover et al.l f2002: 
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iBlecha et al.1 12000). The second epoch data were reduced using 
the standard ESO reduction pipeline for GIRAFFE spectra (ver- 
sion 6.2. a2). Both pipelines apply the standard reduction pro- 
cesses of bias and dark subtraction, and scattered light correc- 
tion. 

The spectra from the 2005 campaign were extracted with the 
optimum method. Those from the 2008 campaign were extracted 
using the standard method. According to the standard extraction, 
for each localized fiber, the spectra are extracted by summing the 
flux of the pixels along the direction perpendicular to the dis- 
persion axis. The optimum extraction method uses the sh ape of 
the fib er profile to weight the flux as a function of noise dHornd 
Il986l) . The use of optimal extraction for the second epoch data 
did not show a sufficient improvement to justify its use. 

For the 2005 spectra, the wavelength solution was calculated 
once at the beginning of the campaign and used for the whole 
data-set. Instrumental drift was followed during the run with si- 
multaneous Th-Ar lamps. For the second epoch data, the wave- 
length calibration was performed using non-simultaneous Th-Ar 
spectra since we did not require as high a precision in radial ve- 
locity as lL08l A sky correction was finally applied to each stellar 
spectrum using standard IRAF routines. 

Combining the two campaigns, the sample of GIRAFFE data 
contains spectra for 1914 different CoRoT targets with a SNR ra- 
tio ranging from 10 to 100. The analysis of such a large amount 
of spectra requires a completely automated and homogeneous 
processing of the whole sample. To render every spectrum com- 
parable to a reference library, we applied the following proce- 
dure. 

The raw images contained a significant number of grazing 
cosmic rays contaminating the dispersion axis. Some spectra 
were polluted by these high values spikes. We corrected the data 
for their presence using a sigma clipping technique adjusted to 
the wavelength domain. This cleaning is a mandatory step to 
avoid numerical problems during the spectroscopic analysis. 

We measured the barycentric radial velocity of the stars us- 
ing a weighted cross-correlation of each s pectrum with a numer- 
ical mask. The latter was constructed by iBaranne et al.l (Q996), 
who identified in the solar spectrum the spectral lines relevant 
to radial velocity measurements. The resulting cross-correlation 
func tion (CCF) is fitted by a Gaussian function. As recalled by 
L08, it provides the barycentric radial velocity (V ra d) and an es- 
timate of th e proje cted rotational velocity (Vsin/). Following 
ISantos et al.l ((2002) methodology, we calibrated the relation be- 
tween the Vsin / and the broadening of the CCF (cr), for the 
FLAMES/GIRAFFE instrument. 

Vsini = Ajcr 2 -crl where a = ^MK (1) 

\ 2V21n2 v ' 

For the slow rotators with typical Vsin / <20 kms -1 , we found 
the coefficient A = 1.8 ± 0.1 by adjusting the slope of Eq. 1 with 
synthetic stellar spectra. The minimal broadening cr in Eq. 1 as 
a function of the (B - V) color index (see Fig. Q] upper panel) 
was fitted by the relation 

cr = 10.907 - 8.361(B - V) + 4.376(B - V) 2 - 0.473(B - V) 3 , (2) 

with the B and V retrieved from Exo-Dat. For faster rotators 
(Vsin / > 20 kms -1 ), the previous relation is no longer valid. 
We convolved a set of 10000 interpolated synthetic spectra 
with different rotational broadening values ranging from 1 .0 to 
80.0 kms -1 . The fit of the Vsin! as a function of the FWHM, 
shown in Fig.[T](lower panel), gives 

Vsin/= -14.22 + 1.04 FWHM + 0.03 FWHM 2 - 0.001 FWHM 3 . (3) 

We applied Eq. 1 to the slow rotators (Vsin/< 20 kms -1 ) and 
Eq. 3 to the faster rotators. We checked that the rotational veloc- 




* * - 
5 r Vs\ni=\.& SQRT((7 z -(7 o 2 ) i 

0.6 0.8 1.0 1.2 1.4 1.6 

B-V 




30 40 50 60 70 80 

FWHM [km.s -1 ] 



Fig. 1. Top: Calibration of the minimal broadening of the 
CCF measured as a function of the (B - V) color index for slow 
rotators (Vsin;< 20 kms -1 ). Bottom: Correlation between the 
V sin i and the FWHM of the CCF for fast rotators. 



ity calculated for stars close to 20 km s 1 gave very close results 
with both methods. 

The derived radial and projected rotational velocities are 
listed in Table [9] A quality flag has been added for V sin i, as 
explained in more detail in Sect. [5] and Table [9] For the SRcOl 
field, no V magnitudes are available in Exo-Dat. We thus ap- 
plied only Eq. 3 to derive the Vsin / value of the stars in the 
SRcOl field, as clearly specified in Table [9] Taking into account 
the limitations of Eq. 3, we preferred to indicate a lower limit of 
40 km s -1 when V sin i > 40 km s -1 . As described in the follow- 
ing section, we applied MATISSE only to the slower rotators in 
the sample since our tests showed that the results are not affected 
by the stellar rotation as long as V sin i is lower than ~ 1 1 km s -1 . 

Each spectrum was then corrected for its V ra d. The same 
analysis with different masks (F0, K5, M4) showed no signifi- 
cant improvements on the velocity precision, 0.2 kms -1 , and 
no important mask effect on the final value. 

We estimated the position of the pseudo-continuum, that is 
the apparent continuum. For that purpose, the spectral lines were 
first removed by an iterative sigma-clipping method. The itera- 
tions were stopped when the dispersion of the remaining pixels 
reached the spectrum noise level. The position of the pseudo- 
continuum was then found by fitting a low order polynomial 
. This first order normalized spectrum was then iteratively im- 
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proved with the results from the MATISSE analysis, as described 
in Sect. 1431 



£ o.io 



center 
anticenter 



-50 50 

Vrad [km.s"'] 




Fig. 2. Kinematics results: Distribution of the barycentric radial 
velocities in the two pointing directions {top) and the projected 
rotational velocities {bottom) for the whole sample. 



following, for high SNR theoretical spectra, we consider the re- 
sults obtained with a direct inversion method. For our observed 
spectra, we noted that the iterative algorithm provided more re- 
liable results. 

The next section describes the grid used for the learning 
phase of MATISSE. We then describe the procedure and estimate 
the errors derived by applying MATISSE on the stellar spectra. 

4.1. The grid of synthetic spectra 

For the training of MATISSE, a grid of theoretical spectra with 
the same spectral resolution and sampling as our observed data 
is required. For the computation of this grid, one has to keep in 
mind that the reliability of the atomic data is crucial when deriv- 
ing accurate parameters from the stellar spectra. We used a list 
of atomic lines derived from the Vienna Atomic Line Database 
(VALD: KupkaetaLlfl99l and lists of molecular lines for the 
species CH, C 2 , CN, OH, MgH, SiH, CaH, FeH, TiO, VO, and 
ZrO and their corresponding isoto pes (kindly provided by B. 
Plez - see iRecio-Blanco et aLlfeOOfl for a complete description). 
We calibrated the atomic line-list with observed spectra of the 
Sun (very high resolution Kuruc z solar spectrum) and Arcturus 
(spectra from lHinkle et al.ll2003l) . 

The parameters used for the spectral synthesis of these stars 
are given in Table [3] an d the solar abundanc es are the same as 
those indicated bv lRecio-Blanco et alj (2006I). For Ar cturus, we 
used the abundances described in Smith et ail (|2000) scaled to 
our solar model. 



Table 3. Atmospheric parameters used for the synthetic spectra 
of the stars used for checking the line-list. 



Star 


T c ff 
(K) 


(dex) 


[M/H] 
(dex) 


[tt/Fe] 
(dex) 


Ref. 


Sun 


5777 


4.44 


0.0 


0.0 


1) 


Arcturus 


4300 


1.7 


-0.6 


0.2 


2) 


Procyon A 


6500 


4.0 


0.0 


0.0 


3) 



References. 1) iGustafsson et al] <2008h . 2) ISmith et al] d2000h . 3) 
Allen de Prieto et al] d2002h 



4. Automatic parametrization of stellar spectra with 
the MATISSE algorithm 

The stellar characterization of our sample includes the determi- 
nation of the following spectroscopic parameters: the effective 
temperature (T e ff), the surface gravity (log g), the overall metal- 
licity ([M/H]), and the c-enhancement ([a/Fe]). 

For that purpose, we used the MATISSE (MAtrix 
Inversion for Spe ct ral Sy nthEsis) algorithm, described in 
Recio-Blanco et al .] (|2006), which is connected to Local 
MultiLinear Regression methods. The stellar parameters are es- 
timated by projections on relevant functions, B{A), derived from 
a multi-linear regression. Since it was initially developed for the 
Gaia-RVS instrument, it is a very efficient means of analyzing ro- 
bustly and automatically large samples of stellar spectra. In the 
present study, we trained MATISSE to the FLAMES instrument 
and in particular to the HR9B setup. 

The algorithm behavior as a function of SNR led us to ap- 
ply an iterative inversion in the computatio n of the B{A) func- 
tions during the learning p hase of MATISSE ( Recio-B lanco et all 
200$ iBiiaoui et al.ll2008l) . in the case of too noisy data. In the 



For the line-list calibration, we first applied the osc illator 
strength modifications described in Gustafs son et alj d2008l) . and 
checked that these modifications improved our fit of the solar 
spectrum, at our working spectral resolution and range. We then 
calibrated the oscillator strengths of about 300 atomic lines at 
our working resolution, while verifying the coherence with the 
high resolution solar spectrum. The^ 2 between the observed so- 
lar spectrum and the theoretical one was improved by a factor 
two. We also calibrated a few lines in the Arcturus spectrum, 
checking that it did not degrade the agreement on the solar spec- 
trum. Finally, we checked that all these modifications also fitted 
the Procy on A spectrum fro m the UVES Paranal Observatory 
Project (Bagnulo et al. 2003 and parameters in Table[3]l. We pre- 
ferred not to modify any oscillator strengths with this star since 
Procyon A parameters and furthermore its abundances are not 
accurately enough known. The modified atomic line data are pre- 
sented in Table [4] No modification was made to the oscillator 
strengths of the molecular data. 

From this calibrated line-list, we computed a grid of theoret- 
ical stellar spectra b ased on a new generatio n of MARCS stellar 
model-atmospheres (Gustafs son et al.l l2008) with the turbospec- 
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Table 4. Shortened list of the modified atomic data. 



Element A (A) Excitation Potential (eV) log gf 



Fel 5129.630 3.943 -1.85 

Fel 5133.681 4.178 0.14 

Fel 5137.382 4.178 -0.4 

Fel 5139.251 2.998 -0.741 

Fel 5139.462 2.94 -0.509 

Fel 5141.739 2.424 -1.964 

Fel 5142.446 4.559 -1.45 

Fel 5142.494 4.301 -0.739 



Notes. The lines not reported here are the values from VALD. The full 
version of this table is available in an electronic table. 



trum code varez & Plezll 19981 and further improvements by 
Plez). The spectra were calculated assuming a plane-parallel ge- 
ometry and % t — 1 kms -1 for stars with log g > 3.5 dex, and a 
spherical geometry, £ f = 2 km s _1 , and a solar mass for star with 
log g < 3.0 dex. The ca lculation of the grid follo wed the same 
process as described by Recio-Bla nco et all (120061) . The final li- 
brary contains 1 1 940 spectra and covers the spectral domain 
from 5 141.70 to 5 347.15 A with a sampling of 0.07 A and a 
resolving power of 25 900. The parameter ranges covered by the 
grid are presented in Table [5] For stars with [M/H] > +0.0 dex, 
we took [a/Fe] = 0.0 dex and for stars with [M/H] < —1.0 dex, 
we chose [a/Fe] = 0.4 dex. Between these two metallicities, a 
linear relation as a function of the metallicity was computed for 
the [or/Fe] value. The [ar/Fe] variations (-0.4, -0.2, 0.0, +0.2, 
+0.4) around this law were then considered for the calculations 
of the spectra. 



Table 5. Stellar atmospheric parameters of the reference grid 
used by MATISSE. 





T e g 




[M/H] 


[or/Fe] 




(K) 


(dex) 


(dex) 


(dex) 


Min 


3000 


1.0 


-3.0 


+0.8 


Max 


8000 


5.5 


+ 1.0 


-0.4 


Steps 


200 or 250 


0.5 


0.25 or 0.50 


+0.2 



To test the validity of our synthetic spectra, we applied the 
whole procedure to the spectra used for the calibration of the 
atomic data. It is mandatory to verify that our procedure is 
self-consistent and that we do not neglect any sources of un- 
certainty. We ensured that the stellar parameters obtained for 
these stars are consistent with the values in Table [3] To test this 
calibration, we also processed the GIRAFFE solar atlar^] in ex- 
actly the same way as our FLAMES observations. The mean pa- 
rameters obtained over the 129 spectra are T e ff = 5740+4 K, 
log g = 4.48+0.01 dex, [M/H] = 0.0+0.004 dex, and [or/Fe] = 
0.02 ± 0.002 dex (the error bars coming from the dispersion). 
These results show the reliability of this calibration and, by ex- 
tension, of the results obtained in this study. The absolute pre- 
cision of our procedure was evaluated by performing a series of 
tests described in the following subsections. 



2 http://www.eso.org/observing/dfo/quality/GIRAFFE/pipeline/ 



4.2. Relative precision of the parameters 

To evaluate the relative precision from one star to the other, 
we applied MATISSE to a large number (10 4 ) of interpolated 
synthetic spectra with randomly chosen parameters. We added 
Gaussian noise to these spectra simulating SNR values repre- 
sentative of our sample: 10, 20, 50, and 100. Figure[3]shows the 
evolution of the error at 70% of the error distribution as a func- 
tion of the SNR. This error is caused only by the dispersion in 
the results since the systematic uncertainty is always negligible. 
These dispersion errors become negligible at SNR > 100. For 
cool metal-rich stars, which exhibit more spectral lines, more 
information about the stellar parameters is available, resulting 
systematically in lower dispersions. The analysis of the cumula- 
tive error distribution (Tableland Fig. [3) showed that most of 
the parameters are recovered with a very high accuracy. 

4.3. Other sources of uncertainty 

We explored the impact on the stellar atmospheric parameters 
derived by MATISSE of a potential error in the radial velocity, 
the normalization, and the effect of the rotational velocity. 

A large uncertainty in the radial velocity will result in a poor 
correction of its effects, creating an error in the determination of 
the various stellar parameters. Several interpolated spectra were 
shifted by different values of the radial velocity to evaluate these 
source of errors. We found that, as long as the error in the V ra d 
remains lower than 1 kms~'(~10% of the resolving power), the 
final uncertainty in the MATISSE parameters is lower than the 
internal precision, presented in Sect. 14.21 Since the mean error 
in the estimated radial velocity for our observed targets is about 
200 ms _1 , we conclude that this source of uncertainty can be 
neglected in our study. 

The rotation of the star broadens the spectral lines, and can 
also introduce an error in the stellar parameters. The MATISSE 
algorithm was trained for non-rotating stars, which represent the 
vast majority of our sample (see Fig. |2J. Using a set of broad- 
ened synthetic spectra, we explored the effect of the stellar ro- 
tation on the photospheric parameters. We found that the pre- 
cision of the parameters is smaller than the internal precision 
when V sin / < 10 kms -1 . Stars with higher rotational velocity 
were disregarded in this study. 

Finally, we investigated the influence of the uncertainty on 
the shape and level of the continuum during the normalization 
process. For this purpose, we calculated a set of theoretical spec- 
tra at a SNR of 20. We modeled a departure from the true contin- 
uum by multiplying each of them by a polynomial of the second 
degree. Each coefficient was fixed at a level of 20% of the er- 
ror in the position, 10% of the error in the slope, and 10% of 
the error in the second order shape of the continuum. Feeding 
these very badly normalized spectra into MATISSE, we found 
that the errors due to the normalization only are AT e ff=175 K, 
Alog g=0.274 dex, A[M/H] =0.388 dex, A[a/Fe] =0.207 dex at 
most (see Sect. 14.4b . This study confirms that the normaliza- 
tion is a matter of prior importance and must be carefully taken 
into account. We noted that the effect is more important for high 
SNR spectra for which normalization is a critical issue. To over- 
come this issue, we implemented supplementary iterations for 
the placement of the continuum. Each observed spectrum was 
divided by the synthetic spectrum calculated with the first es- 
timate of the photospheric parameters obtained with MATISSE. 
The residual of the division was then cleaned from remaining 
spectral lines and fitted by a low order polynomial. The observed 
spectrum was divided by this new continuum. A new estimate of 
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Table 6. Maximum absolute value of the internal error for various proportions of the interpolated spectra at SNR=10, SNR=20 and 
SNR=50 







SNR=50 






SNR=20 






SNR=10 






70% 


80% 


90% 


70% 


80% 


90% 


70% 


80% 


90% 


T eff (K) 


21 


30 


50 


50 


65 


98 


76 


105 


175 


log g (dex) 


0.03 


0.04 


0.06 


0.08 


0.10 


0.14 


0.13 


0.17 


0.24 


[M/H] (dex) 


0.02 


0.03 


0.04 


0.05 


0.07 


0.09 


0.08 


0.10 


0.14 


[ar/Fe] (dex) 


0.01 


0.01 


0.02 


0.02 


0.03 


0.04 


0.04 


0.06 


0.08 



the photospheric parameters was made using this re normalized 
spectrum. A few iterations were sometimes necessary before 
reaching a stable solution for the atmospheric parameters. This 
iterative normalization ensures the minimization of this source 
of uncertainty that can cause large errors in the estimation of the 
stellar parameters. 



4.4. Comparison with libraries of stellar atmospheric 
parameters 

To validate the stellar parameters derived with MATISSE, we 
retrieved spectra from various libraries: th e Elodie 3.1 archive 
dPrugniel et al I 120071). the S 4 N study dAllende Prieto et all 
120041) . and fhe lSantos etail ((2009) work. We adapted these spec- 
tra to our spectral domain and resolution. 

The Elodie library contains a quality flag for the T e ff, log g, 
and [Fe/H] taken from the literature. To ensure that the de- 
termination of the spectroscopic parameters are of good ac- 
curacy, we selected only the stars with high quality criteria 
for the T e ff and the metallicity. These stars presented a poor 
log g quality criterion, so we decided not to include it in 
our sample. We also selected 90 stars from the S 4 N study, 
the remaining 29 being either fast rotators or stars for which 
the parameters were not determined in this study. The sur- 
face gravity provided in the S 4 N study was not derived spec- 
troscopically. The spectroscopic method used by these au- 
thors showed a too high dis crepancy with literature parame- 
ters dAllende Prieto et al. 200j). He nce, this parameter was de- 
rived by the algorithm described by All ende Prieto et al. I J200l 
from isochrones tracks and parallaxes from the HIPPARCOS 
catalogue. For this sample, we therefore compare spectroscopic 



Table 7. Uncertainty at 70% of the error distribution on the mea- 
surement of stellar parameters with spectra from the literature. 



gravities from MATISS E to evolutionary ones. In the literature, 
and in particular in Allende Priet o et all (120041) . the metallicity 
is inferred from the iron content ([Fe/H]). If the star appears 
not to have peculiar abundances, this value should be equal to 
the overall metallicity. For the [a/Fe] comparison, we used the 
abundances fo r the neutral magnesium, c alcium, and silicium, 
determined by lAllende Prieto et al" (2004), since these lines are 
presen t in our spectral range. Finally, the study of Santos et al. 
(2009) was interesting because it includes giant stars. We se- 
lected only these stars so as to probe the efficiency of MATISSE 
for giant targets, uncovered by the two previous studies. 

In Fig. |4] we compare MATISSE spectroscopic parameters 
as a function of the literature ones. The parameters derived by 
MATISSE are available in electronic Tables [10] \U\ and [13] We 
find very good agreement over the four atmospheric parameters: 
70% of these stars have difference lower than ~85 K in T e ff, 
-0.2 dex in log g, -0.15 dex in [M/H], and -0.1 dex in [a/Fe] 
(see Table[7]l. This shows that for these high SNR observed spec- 
tra, the total error on the estimation of the parameters is at least 
four times higher than for interpolated synthetic spectra (inter- 
nal error negligible for SNR > 100 — see Sect. 14.21 and Fig. 



Nb. Spec. 


T e ff 


l°g£ 


[M/H] 


[a/Fe] 


Ref 




(K) 


(dex) 


(dex) 


(dex) 




118 


61 




0.11 




1) 


90 


84 


0.202 


0.149 


0.079 


2) 


39 


75 


0.184 


0.129 




3) 



References. l)lPrugniel et al.l d2007l) . 2) 1 Allende Prieto et aljj2004b . 3) 
ISantos et al] <2009l) (see al"soFig.l4b 



0. This illustrates that the total error in the estimation of the 
parameters is caused not only by the method, which exhibits an 
external source of uncertainty. 

We took advantage of multiple observations in our sample to 
evaluate the real total error on the estimation of every spectro- 
scopic parameter. About fifty stars in the center direction were 
observed twice and analyzed separately by the entire pipeline. 
The SNRs of these data range between 5 and 60. The error at 
70% of the distribution of the difference between the two deter- 
minations provides another estimate of the uncertainty. We com- 
bined this uncertainty with the internal error and the external er- 
ror (quadratic sum) to estimate the total errors of <x Xdr = 140 K, 
o"iog g = 0.27 dex, ct[ M /h] = 0.19 dex, and cr^Fe] = 0.09 dex. 
We recall that for these samples, the relative precision from one 
star to another is much lower and given in Table [6] 

5. Results 

Table [8] reports a brief summary of the numbers of C0R0T stars 
presented in our study. It lists the number of observed targets, the 
number of stars for which we derived the physical parameters 
with MATISSE, the number of detected double-lined spectro- 
scopic binaries (SB2), and the number of stars with estimated ra- 
dial and projected rotational velocities. For the targets observed 
in the LRaOl, we used the radial velocities as reported in the 
study performed by L08. Most of the targets presenting no CCF 
in the center direction were missed due to the poor quality of the 
input astrometry in the SRcOl field. 

The barycentric radial velocities, presented in Table [9] were 
measured for 77% of the stars in the center direction and 85% 
stars in the anticenter direction. We only measured precise radial 
velocities for stars presenting a CCF with a reasonable contrast 
and not dominated by strong rotational broadening. The V sin/ 
value was measured for 1 604 stars in our sample (Table [9), in- 
cluding objects for which V ra d had not been well determinecQ. 
As mentioned in Sect. [3] the V sin i estimations have a quality 
flag, reported as an exponent of the value in Table [9] Its value 
is set to be 1 if the V sin i estimate is not affected by noise or a 



3 In this case, the Gaussian can still be fitted on the wings of the CCF. 
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Fig. 5. Hertzsprung-Russell diagram with the T e g and log g de- 
rived by MATISSE. The red crosses are stars from the LRcOl and 
the blue triangles from the LRaOl field. 



Fig. 6. Distribution of the overall metallicity in the two pointing 
directions for giant (left) and dwarf (right) stars, normalized to 
the number of stars in each field. 



second component in the CCF, 2 if the target is a SB2 and the 
V sin i is related to the main component of the CCF, and 3 if the 
contrast and shape of the CCF are insufficient to assure a proper 
estimate of its parameters. As shown in Fig. [2] most of the stars 
are slow rotators (V sinj < 10 k ms 1 for 75%) as expected for 
late-type field stars (Grav 2008). In both samples, we excluded 
from the MATISSE analysis the stars exhibiting a CCF with a 
FWHM greater than 20 kms -1 , SB2, and stars for which no ac- 
curate radial velocity was measured because of the bad quality 
of the CCF. 

The MATISSE spectral analysis, described in the previous 
section, was carried out on a total of 1 227 targets that is 65% 
of the initial number of stars. The results of this analysis are 
presented in Table [9] Combining the T e ff and log g values, we 
compiled the Hertzsprung-Russell diagram presented in Fig. [5] 
We observed main-sequence solar-type and cooler giant stars in 
each of the CoRo T/Exoplanet fields. 

We used the log g derived by MATISSE to differentiate gi- 
ants, log g < 3 .1 dex, and dwarfs, log g > 3 .8 dex. This criterion 
is based on the Stra izvs & KuriheneT(fl9 8 1 ) tables. The metallic- 
ity distribution for the giants and dwarfs are presented in Fig. [6] 
It is similar in the two pointing directions for the dwarfs. For the 
giants, the maxima of the metallicity distributions are apart one 
from the other of a few tenth of dex. This could be explained by 
the r adial metallicity grad ient generally observed in the Galactic 
disk (Ped icelli et al.ll2009l and references therein). We compare 
in Fig. [7] the dwarf-giant dichotomous separation based on the 
log g derived by MATISSE with the color-magnitude diagram r' 



Table 8. Number of targets according to the kinematics analysis. 



Number of 


LRcOl 


SRcOl 


LRaOl 


Total 


Stars observed 


689 


484 


741 


1914 


Stars analyzed by MATISSE 


555 


215 


457 


1227 


SB2 detected 


17 


7 


17 


41 


V rad determinations 


630 


270 


634 


1534 


V sin i estimates 


658 


268+19 


659 


1604 



Notes. For the slow rotators in the SRcOl field, we could not apply Eq. 1 
from Sect. [3] since no V magnitude is available in Exo-Dat. We applied 
Eq. 3 for these 268 targets. 



versus (r'-J). According to this color-magnitude diagram, the 
stars with an intermediate log g value ([3.1-3.8]) could be con- 
sidered as giants. We note that the dwarf-giant dichotomy used 
for this first order classification does not include the uncertainty 
in the log g. In the following, we include these stars in the giant 
sample. This figure also shows how the spectroscopic classifica- 
tion easily allows us to distinguish the dwarf and giant popula- 
tions, this issue being one of the most constraining limitation of 
the photometric classification, especially for faint targets. 

We find a good agreement between these simple photomet- 
ric and spectroscopic classifications for the majority of our sam- 
ple. However, a few stars are misclassified in both pointing di- 
rections. In the center field, the outliers could correspond to 
reddened dwarfs. In the anticenter direction, the obvious mis- 
classification is due to the limits of the grid of the theoretical 
spectra (see Sect. 14. U . This object is a late-A metal poor star 
([M/HN -0.7 dex) that exhibits very few lines in the spectral 
range, resulting in a very difficult classification. We found that 
~ 1% of the stars are misclassified. 

The photometric classification that is available in Exo-Dat 
consists of two steps. First, a raw luminosity class was de- 
rived from color-magnitude diagrams similar to Fig. [7] Then 
the final estimates of luminosity classes and spectral typ es were 
based on spectral energy distribution (SED) fitting with Pickles 
( 1998) templates. For a full description, we refer to lDeleuil et al. 
(2009). No preliminary classification was available for stars in 
the SRcOl field, hence the center simply stands for the LRcOl 
field. In Fig. [8] we compare the results of the spectral anal- 
ysis performed by MATISSE with the photometric classifica- 
tion available in Exo-Dat. The conversion from T e ff-log g to 
spectral type and luminos ity class was performed using the 
Straizvs & Kurilieni d!981l) tables. Figure [8] illustrates the rea- 
sonable agreement found for the spectral types, whereas the lu- 
minosity classes histograms show clearly more dwarfs and less 
giants than estimated by the photometric classification. A dis- 
crepancy was expected between the photometric and spectro- 
scopic luminosity classes since the initial photometric estimation 
of the luminosity class consisted of a simple cut-off in a color- 
magnitude diagram, applied to prevent the omission of dwarfs 
during the CoRoT target selection phase. This cut-off could re- 
sult in a misidentification of dwarfs and giants. For the faint tar- 
gets, this becomes a more critical issue since the two populations 
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are mixed (Fig. |7J. We indeed found a closer agreeement for the 
brightest targets in our sample. 

6. Discussion 

The samples of CoRoT targets for which we determined atmo- 
spheric parameters are not representative of the whole fields ob- 
served by CoRoT. As described in Sect. [2] we used selection 
criteria for the multi-fiber observations. As a result, our analy- 
sis was focussed primarily on the slowly-rotating population of 
the FGK stars. To more accurately characterize the stellar pop- 
ulation in the exoplanet fields, we took into account the vari- 
ous biases introduced by our selection. We de biased our sam- 
ples using density plots of the color-magnitude diagrams pre- 
sented in Fig. |7] The limitations of this study are constrained by 
the ranges in color and magnitude covered by our observations 
(blue boxes in Fig. [7j. For every bin of 0.5 in r' and 0.5 in (r'- 
J), we randomly selected stars in our observations so as to fill 
in the two-dimensional distribution in the same proportions as 
the Exo-Dat one. The number of stars in the two long run fields 
was large enough to efficiently apply this procedure. This was 
not the case for the SRcOl field to which we did not apply this 
de-biasing procedure, so the center field simply stands for the 
LRcOl field. To validate this procedure, w e took advantage of 
the Besangon Galactic model (Ro bin et al.ll20 03). We compared 
each de-biased field with a simulation taking into account the 
limits in magnitude and assuming the standard extinction law. 
We found reasonable agreement for the distributions of T e g and 

To compare the two de biased samples, we selected in LRcOl 
only stars with r 1 < 15 mag since this is the magnitude limit of 
our observations in the LRaOl field. For the two fields, we in- 
cluded the stars with an intermediate log g in the giant sample 
since they are too few to be analyzed separately and they be- 
long mainly to the giant branch in the color-magnitude diagram 
(Fig.pJ. 

Figure [9] shows the resulting distributions of the T e ff, [M/H], 
and [c/Fe] derived by MATISSE for the two fields normalized 
to the entire de biased sample. We used the log g values from 
MATISSE to separate giants from dwarfs as described in the pre- 
vious section. The dwarf content is roughly similar in the two 
directions over these three parameters. The majority of these 
stars are early G or late F stars with solar metallicity. The LRcOl 
field clearly contains mainly late K giant stars. The other field 
presents a smoother distribution around K giants, which are 
slightly hotter. The histograms of [M/H] and [a/Fe] for the gi- 
ant stars ill ustrate the metallicity gradient generally observed in 
our galaxy dPedicelli et al.ll2009h . with the center field contain- 
ing more higher [M/H] giant stars than the anticenter one. For 
targets with r' < 15 mag, the two fields can be compared: the 
LRcOl field is composed of about 41% giant and 59% dwarf 
stars, whereas the LRaOl field is composed of about 26% giant 
and 74% dwarf stars. Adding the objects with r' > 15 spectro- 
scopically observed in the center direction, we estimated the to- 
tal amount of dwarfs to be 72% of the stars observed by CoRoT 
in this field. 

To estimate the detection rate of exoplanets in the LRcOl 
and LRaOl fields, we applied the probability law from 
lUdrv & Santol (l2007h {V = 0.03 x io 204x[M/H] for [M/H] > 
0.0 dex and V = 3% for [M/H] < 0.0 dex) to the FGK dwarf 
population with r' < 15. According only to this metallicity cri- 
terion, we expected 3.6% of planets in LRcOl field, and 3.7% in 
the LRaOl field. The geometrical probability that a transit occurs 
is V — R*/a, where a is the semi-major axis of the planetary orbit 



and R„ the host sta r's radius. We use d the distribution of semi- 
major axis from the ISchneiderl (I 1 9951) website and integrated the 
transit probability over this distribution. Combining these prob- 
abilities and applying them to the two selected fields, we found 
that the number of expected transiting planets detection is 8+7 
and 9±7 for the LRcOl and LRaOl fields, respectively, with the 
error bars simply resulting from the propagation of the error in 
[M/H] (see Sect. [44}. 

This results is consistent with those derived so far by the 
CoRo r/Exoplanet consortium, i.e. four transiting planets in the 
LRcOl field and three in the LRaOl field. These estimates must 
be interpreted with caution because the metallicity probability 
relation assumed above was derived from planetary detections 
made using radial velocity techniques whereas the CoRoT space 
mission bases its detection on planetary transits. 



7. Summary 

Using the MATISSE algorithm, we have developed a fully au- 
tomated pipeline to perform a homogeneous spectral analysis of 
FLAMES/GIRAFFE spectra covering the Mg I b line wavelength 
range with an intermediate resolving power (HR9B setup). We 
have estimated the atmospheric parameters (T e ff, log g), metal- 
licity indices ([M/H], [a/Fe]), the barycentric radial velocity 
( V ra dX an d the projected rotational velocity (V sin i) of more than 
a thousand of stars located in three of the CoRo T/Exoplanet 
fields towards the Galactic center and anticenter: the LRaOl, 
LRcOl, and SRcOl fields. 

Comparing these results with the photometric classification 
available in Exo-Dat, we have found reasonable agreement for 
the spectral types and a small discrepancy between the luminos- 
ity classes. This is due to the difficulty in deriving a luminosity 
class based only on the broad-band photometry available in Exo- 
Dat, since this consists of a simple cut-off in the color-magnitude 
diagrams (see Fig. O. This limit becomes very uncertain for very 
faint targets (r' > 14.5). 

From the derived stellar atmospheric parameters, we have 
compiled the first un-biased reference sample for studying the 
planet-metallicity relation in the CoRo T/exoplanet fields. We 
found that the number of planets discovered so far by CoRoT 
is in agreement with the probability law predicting the number 
of planets as a function of the stellar metallicity. 

In the near future, we will easily be able to implement these 
efficient procedures to any other CoRoT field. FLAMES obser- 
vations are currently being performed in other CoRoT pointing 
directions: the LRcOl and LRc03. The spectral analysis of 1 000 
stars with MATISSE can be performed within a few minutes us- 
ing an average computer. This would increase the statistics of 
available data and help us to understand the typical stellar pop- 
ulation observed by CoRoT. Comparing the parameters of the 
planet-hosting stars with the more general stellar field they be- 
long to, will provide additional constraints on planetary system 
formation scenarios. 

Using the atmospheric parameters with stellar evolution 
models will lead to the determination of physical parameters for 
a large number of CoRoT targets. By adding proper motions to 
the parameters derived in this work, we will be able to derive 
kinematical information about these sample stars and help in ex- 
ploring the age-metallicity and age-kinematics relationships in 
these fields. Combining this key information with the richness of 
data provided by the light curves would provide unprecedented 
insight into these relationships. This will be explored in a forth- 
coming paper. 
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The spectroscopic parameters, the kinematics information 
derived, and the spectra used in this study will be made avail- 
able to the community through Exo-Dal^ . 
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Fig. 3. Internal error at 70% of the error distribution as a function of the SNR for cool dwarf metal-rich stars (T e ff < 5500 K, 
log g> 4.0 dex, [M/H] > 0.0 dex, green diamonds), cool giants with intermediate metallicity (T e ff< 5500 K, log g< 3.0 dex, 
0.0 > [M/H] > -1.5 dex red stars), hot subgiants metal poor (T eff > 5500 K, 4.0 >log g> 3.5 dex, [M/H] < -0.75 dex, black 
crosses), and the whole sample of interpolated spectra (blue triangles). 

Table 9. Derived physical parameters for the stars observed in the LRaOl, LRcOl, and SRcOl CoRoT fields. The error bars for the 
atmospheric parameters are only internal errors coming from Fig. [3] The V sin i error is of the order of 10% of its value. The quality 
flag (exponent of the value) was set to 1 if the V sin i estimate is not affected by noise or a second component in the CCF; 2 the 
target is a SB2 and the V sin i is related to the main component of the CCF; 3 the contrast and shape of the CCF are insufficient to 
assure a proper estimate of its parameters. The SNR of the FLAMES/GIRAFFE spectrum is listed in the last column of the table. 
We remind the reader that the coordinates, observing dates and magnitudes are publicly available through Exo-Dat. 



CoRoT ID 


V sin i 


v rad 




T e ff 




logs 


C"log g 


[M/H] 


O" [M/H] 


[o/Fe] 


C [a/Ft] 


SNR 




(km s" 1 ) 


(kms- 1 ) 


(kms" 1 ) 


(K) 


(K) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 




100603128 


3.8 1 


-7.83 


0.187 


5350 


79 


4.01 


0.13 


-0.24 


0.08 


0.06 


0.04 


9 


100583300 


17.8 1 


-61.10 


0.163 


















23 


100567221 


1.4 1 


-25.66 


0.075 


5889 


48 


4.29 


0.08 


-0.04 


0.05 


0.02 


0.02 


23 


100565715 


5.8' 


92.85 


0.246 


6445 


71 


4.64 


0.12 


-0.25 


0.07 


0.16 


0.03 


12 


100577769 


8.2 1 


-27.26 


0.272 


5117 


76 


3.02 


0.13 


-1.02 


0.08 


0.24 


0.04 


10 


100588558 


11. 3 1 


39.90 


0.299 


















13 


100576988 


4.8 1 


-26.10 


0.118 


5963 


60 


4.41 


0.1 


0.0 


0.06 


-0.02 


0.03 


16 


100636016 


2.7 1 


22.55 


0.322 


5673 


80 


4.29 


0.13 


-0.88 


0.08 


0.24 


0.04 


8 


100637229 


19.4 1 


-12.39 


0.224 


















28 


100604545 


10.7 1 


0.37 


0.231 


















10 


100634976 


3.9 1 


-37.27 


0.221 


6100 


73 


4.07 


0.12 


-0.35 


0.08 


0.07 


0.04 


11 
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Fig. 4. Comparison of th e stellar atmospheric param eters derived by MATISSE and those of the literatu re. Red crosses : Prueniel et alJ 
(120071) : gr een diamond : lAllende Prieto et alJ d2004l) : blue circles: UVES spectra & parameters from lPrugniel et al.l (120071) : yellow 
triangles : ISantos et al] d2009l) . 



Table 10. Derived parameters for the S 4 N sample from lAllende Prieto et alJ ([2004) (complete table available at the CDS). The first 
line is the Sun. 



Hipp. Num. 


T e fl 
(K) 


(de)] 


[Fe/H] 
(dex) 


[o/Fe] 
(dex) 


T eS MATISSE 
(K) 


log g MATISSE 
(dex) 


[M/H] MATISSE 
(dex) 


[a/Fe]MATISSE 
(dex) 




171 


5777. 
5361. 


4.437 
4.610 


0.000 
-0.770 


0.000 
0.343 


5724 
5376 


4.227 
4.359 


0.036 
-0.631 


0.029 
0.429 



Table 11. Derived parameters for the Elodie3.1 sample from Prug niel et al.l d2007l) (complete table available at the CDS). 



HD identifier 


T c ff 
(K) 


(dex) 


[Fe/H] 
(dex) 


T^MATISSE 
(K) 


log gMATISSE 
(dex) 


[M/H] MATISSE 
(dex) 


[a/Fe]MATISSE 
(dex) 


HD000400 
HD001835 


6146. 
5777. 


4.090 
4.450 


-0.280 
0.170 


6113 
5883 


4.397 
4.804 


-0.185 
0.197 


-0.075 
-0.068 
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Fig. 7. Color magnitude diagrams for the giants (red stars) and the dwarfs (green crosses), as well as the stars with an intermediate 
log g (blue triangles). All the stars in the Exo-Dat (black dots) are over plotted in the diagram. The dichotomy giant-dwarf comes 
from the MATISSE parameters. The blue boxes correspond to the limits in color and magnitude for the de-biasing of our sample. 
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Fig. 8. Comparison of the spectral types and luminosity classes derived with MATISSE and with the Exo-Dat photometric classifi- 
cation criteria. 
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Fig. 9. De-biased distributions of the effective temperature, overall metallicity and a-enhancement in the LRcOI (solid red) and 
LRaOl (dashed blue) fields. The top and lower panels show the distributions for the dwarf and giant stars, respectively. The distri- 
butions have been normalized to the number of stars in each field. The percentage given in the [a/Fe] histogram represents the rate 
of dwarfs (upper panel) and giants (lower panel) in the two fields. 

Table 12. Derived parameters for the UVES (POP. lBagnulo et al.ll2003l) sample and literature from lPrugniel et all d2007l) . 



HD identifier 


T e ir 


logs 


[Fe/H] 


T eS MATISSE 


log gMATISSE 


[M/R]MATISSE 


[a/Fe]MATISSE 


(K) 


(dex) 


(dex) 


(K) 


(dex) 


(dex) 


(dex) 


HD022049 


5089.0 


4.55 


-0.10 


5054 


4.542 


0.010 


0.043 


HD022484 


5989.0 


4.10 


-0.07 


5916 


3.891 


-0.040 


0.024 


HD030562 


5860.0 


4.03 


0.17 


5854 


4.053 


0.271 


0.009 


HD061421 


6562.0 


4.08 


-0.02 


6366 


3.990 


0.067 


-0.085 


HD076932 


5850.0 


4.05 


-0.93 


5600 


3.593 


-1.052 


0.391 


HD1 15383 


5979.0 


4.14 


0.08 


6063 


4.470 


0.229 


-0.047 


HD128167 


6782.0 


4.32 


-0.40 


6658 


4.543 


-0.214 


-0.022 



Table 13. Derived parameters for the Santo s~et~a"fl d2009l) sample (complete table available at the CDS). 



T cff logg [Fe/H] T^MATISSE log gMATISSE [M/H] MATISSE [a/Fe]MATISSE 
(K) (dex) (dex) (K) (dex) (dex) (dex) 



Identifier 



IC2714NollO 
IC2714No87 



5017. 
5029. 



2.850 0.010 
2.620 -0.060 



5028 
5009 



2.911 
2.889 



-0.033 
-0.181 



0.047 
-0.024 
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